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Abstract 

A generalization of Hurwitz stable polynomials to real rational functions is considered. We estab- 
lishe an analogue of the Hurwitz stability criterion for rational functions and introduce a new type of 
determinants that can be treated as a generalization of the Hurwitz determinants. 

Introduction 

It is well known that the problem of stability of a linear difference or differential system with constant 
coefficients reduces to the question of locating the zeroes of its characteristic polynomial in the left half- 
plane of the complex plane. One of the most famous results from stability theory is the Hurwitz theorem, 
which expresses stability of a real polynomial in terms of its coefficients [5] El E] (see also [3IZ])- Namely, 
the Hurwitz theorem states that a real polynomial has all its zeroes in the open left half-plane if and only 
if some determinants constructed with the coefficients of the polynomial are positive (see Theorem II. 5|) . 
Those determinants are now called the Hurwitz determinants due to Adolf Hurwitz who introduced them 

In the present work we investigate a class of real rational functions satisfying the Hurwitz conditions: 
the Hurwitz determinants constructed with the coefficients of the Laurent series at oo of a real rational 
function are positive up to the order n, where n is the sum of degrees of the numerator and denominator 
of the rational function. This is a generalization of the class of real Hurwitz stable polynomials. It turns 
out that such class of rational functions is characterized by location of poles and zeroes: all zeroes lie in 
the open left half-plane of the complex plane while all poles lie in the open right half-plane (Theorem l2.2p . 
Finally, we express the Hurwitz determinants of real rational functions in terms of the coefficients of the 
numerator and denominator (Lemma l2.4l) . Thereby we introduce a new type of determinants and describe 
the class of real rational functions satisfying Hurwitz conditions in terms of coefficients of their numerator 
and denominator (Theorem I2.5[) . As well as in case of polynomials this class of rational functions is 
characterized by positivity of those determinants. Some simple properties of the new type of determinants 
that can be treated as a generalization of the Hurwitz determinants are considered. 
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1 Hurwitz polynomials 

Consider a real polynomial 

p{z) d = a z n + a 1 z n - 1 H h 

Throughout the paper we use the following notation 



01, . 



,a n £ 



a > 0. 



(1.1) 



(1.2) 



where n = degp, and [p] denotes the largest integer not exceeding p. 
The polynomial p can always be represented as follows 

p( z ) =Po(z 2 ) + zp l {z 2 ), 

where po and p\ are the even and odd parts of the polynomial, respectively. Introduce the 
function: 

dg p_i(u) 

pom 

Definition 1.1. We call $ the function associated with the polynomial p. 
Associate with the polynomial p the following determinants: 
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(1.3) 



Ai(p) 



where we set a; = for i > n. 
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Definition 1.2. The determinants Aj(p), j = 1, . . . ,n, are called the Hurwitz determinants or the Hurwitz 
minors of the polynomial p. 



Suppose that degpo ^ degpi and expand the function $ into its Laurent series at oo: 



\ Vim 
$ («) = — r4 =s_i 
Po(«) 
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(1.5) 



where s_i ^ if degpo — degpi, and s_i = if degpo > degpi. 
For a given inhnite sequence (sj)°^ , consider the determinants 
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.7 = 1,2,3,. 



(1.6) 



These determinants are referred to as the Hankel minors or Hankel determinants. 

Together with the determinants (|1.6[) we consider one more sequence of Hankel determinants. 



D 3 m dg 
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Sj+1 



S2j-1 



J = 1,2,3,. 



(1.7) 



It is very well known (5j [3] (see also [7]) that there are relations between the determinants Dj($), 
Dj(<&) and the Hurwitz minors Aj(p): 
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1) If n = 21, then 



2) If 7i = 2l + 1, then 



A 2j - 1 (p) = a 2 j - 1 D j (<i>), 

j = l,2,. ..,1- (1.6 

A 2j -(p) = (-1)%^($), 



A 2j (p) = Q^) 



j = 0,l,...,Z; (1.9) 



where A)($) = 1- 
Here I is dehned in (|1.2[) . 

It is also well known [SJ[3] (see also [3]) that the number of poles of the function <£> equals the order of 
the last non-zero minor Dj(&). Since the function $ has at most I poles, we have 

=£,(#) = (), j>l. (1.10) 

Thus, in the sequel, we deal only with the determinants ZXj($), Dj($>) of order at most i. 

Definition 1.3. The polynomial p defined in (|l.ll) is called Hurwitz or Hurwitz stable if all its zeroes lie 
in the open left half-plane of the complex plane. 

The following criterion of Hurwitz stability of a real polynomial was (implicitly) established in [5] (see 

also Ganid). 

Theorem 1.4. Let a real polynomial p be defined by (|1.1[) . The following conditions are equivalent: 

1 ) the polynomial p is Hurwitz stable: 

2) the following hold 

s_i > for n = 2l+l, 

>0, j = I, (1.11) 

(-1)^($)>0, j = l,..., I, 



where I 



~n~ 
-2- ' 

This theorem together with fornlae (jl.8|l - (|1.9|l imply the following theorem which is very well known 
as the Hurwitz criterion of polynomial stability. 

Theorem 1.5 (Hurwutz [5]). A real polynomial p of degree n as in (jl.lj) is Hurwitz stable if and only if 
all Hurwitz determinants Aj (p) are positive: 

Ai(p) > 0, A 2 (p) > 0, . . . , A n (p) > 0; (1.12) 

2 Hurwitz rational function 

Consider a rational function 

def h(z) = t ^ r _ m + ^r-m-1 + + _ _ _ (2 ^ 

9{z) 



where h and g are real polynomials 



h(z) d = b z r + 6iz r_1 + . . . + 6 r _iz + b r , b , &i, . . . , b r G K, b > 0; 

g( Z ) = C Q Z m + C X Z m ~ l + . . . + Cm-iZ + C m , C , Cy, . . . , C m € R, C > 0. 



(2.2) 



Here r, m € N1J{0}, r + m > 0. Assume that the polynomials p and q are coprime. The number n = r + m 
is called the order of the function R. 
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Definition 2.1. The real rational function R is called Hurwitz rational function if both polynomials h(z) 
and g(—z) are Hurwitz stable- 
It turns out that one can establish a criterion for real rational functions to be Hurwitz similar to the 
Hurwitz criterion of polynomial stability 

Theorem 2.2. A real rational function R{z) of the form (|2.ip - (l2.2l) is Hurwitz if and only if the leading 
principal minors Ak(R) of the infinite Hurwitz matrix 

\ 



U{R) = f 
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(2.3) 



are positive up to the order n (inclusive) : 

Ai(fl) > 0, A 2 (R) > 0, . . . , A n {R) > 0. 
Proof. Consider the following auxiliary polynomial 



def 



P(z) = (—l) m h(z)g(—z) = aoZ n + axZ n + . . . + a n -\z + a n , ao = foo c o > 0. 

By definition, the function R is Hurwitz if and only if the polynomial P is Hurwitz stable. 
Consider the rational function $ associated with the polynomial P: 
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(2.4) 



(2.5) 



where the polynomials Po(u) P\{u) are the even and odd parts of the polynomial P, respectively. Since 
2P (z 2 ) = P(z) + P(-z) and 2zP (z 2 ) = P(z) - P(-z), we have 



z$(z 2 ) 



Hence we obtain 



where 



P(z) - P(-z) = h(z)g(-z) - h{-z)g(z) = R(z) - R(-z) 
P(z) + P(-z) h{z)g(-z) + h{-z)g[z) R(z) + R{-z)' 

Hu) _ Rl (u) _ Pl (u) 



Ro(z 2 ) 
Ri(z 2 ) 

are the even and odd parts of the function R, respectively. Thus, for n = 21 + 1, we have 
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(2.6) 



$(m) 

and for n = 21, we have 

$(u) = 



t + t 2 u 1 + t±u 2 + 
h + tzu- 1 + t 5 u~ 2 + 



hu- 1 + t 3 U- 2 + tgU- 3 + 
t + t 2 U~ 1 + t A U^ 2 + . . 
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We recall that n = r + m is the order of the rational function R. Now note that the formulae (|1.8p - (|1.9p 
are formal, so following verbatim the proof of the formulae (|1.8p - (|1.9p (see, for instance, [TjH]) one can 
establish the following relationships 



1 Or one of them is Hurwitz stable while the second one is a constant. 
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1) for n = 21, 

2) for n = 21 + 1, 



A 2j _ 1 (R) = t 2 j - 1 D j (<f>), A 2j {R) = t$D j {<S>) 5=1,2,...,/, 



2J ! 



to 



2j 



A 2i (i?)= (^-] £>,($), 5 = 1,2,...,/, 



A 2i+1 (i?) = (-l)i 



to 



2j+l 



5 = 0,1,...,/. 



where the Hurwitz minors Aj(R) are defined as follows 
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(2.7) 



(2.8) 



(2.9) 



By Theorem ll.41 the polynomial P is Hurwitz stable if and only if the minors Dj(§) and Dj($>) satisfy 
the inequalities (ll.lip . Now since the polynomial P is Hurwitz stable if and only if R is a Hurwitz rational 
function, the formulas (|2.7p ~ (|2.8p and Theorem 11.41 imply the assertion of the theorem. □ 



Note that the leading principal minors Aj (R) of order more than n equal zero regardless whether the 
function R being Hurwitz or not: 

A n+1 (R) = A n+2 {R) - A n+3 (R) = • • • = 0. (2.10) 

It follows from (11.10[) and from the folrmulas (|2.7j) - p.8j) which are obviously valid for j > I. 



Remark 2.3. It is easy to see that ao = &o c o an d to = — ■ So the formulas (| 1 .8|) — (| 1 .9[) and the for- 

co 

mulas (|2.7p ~ (|2.8p imply the equalities 



A j (P)=cl j A j (R), 
This equalities verify (|2.10[) . since Aj(P) = for j > n. 



J 



1,2,. 



(2.11) 



Recall now that one of necessary conditions for the polynomial p defined in (II. lj) to be Hurwitz stable 
is the positivity of its coefficients a,j > 0, j = 0, 1, . . . , n. This is the so-called Stodola theorem |3]. It 
turns out that positivity of the coefficients tj in (|2.1[) is not necessary condition for the Hurwitzness of the 

function R. Indeed, if R = — is Hurwitz, then the polynomial g has all zeroes in the open right half-plane, 

9 

so it has coefficients of different signs. Therefore, it is easy to find polynomials h and g such that R has 
both negative and positive coefficients. For example, the following function is Hurwitz, but its Laurent 
series at infinity has positive, negative and zero coefficients: 



F{z) 



z z + z + 1 



t + t lZ 1 +t 2 z 2 + t 3 z 3 + .... 



where 



and 



h 3 -2 =t 3j -i = (-i) y - x) 2, t 3j 







for 



J 



1,2,. 



Finally, we find a connection between the coefficients of the polynomials h and g and the Hurwitz 
determinants Aj(R). This gives us a criterion of Hurwitzness of a real rational function in terms of the 
coefficients of its numerator and denominator. 
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Let again the function R be defined by (|2.f p ~ (|2.2p . Introduce the following determinants of order 2j, 
i = 1,2,...: 



Q 2j (h,g) 
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(2.12) 



Here we set bk '■= for k > r, and Cj := for j > m. 
Lemma 2.4. For any j — 1,2, ... , 



n 2j (/i, ff ) = c^A i ( J R), 
where Aj (i?) are i/ie Hurwitz determinants of the function R defined in (|2.9 
Proof. First interchange the rows of the determinant (|2.12p to obtain 



(2.13) 
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(2.14) 



This does not change the sign of tt 2 j(h,g). In fact, lower the jth and (j + l)st rowf@ to their initial 
positions in (|2.12p . This will require an even number of transpositions. The next pair of rows will then 
meet, the lowering operation will be applied to them, and so on (J — 1 times) until we obtain the initial 
determinant (|2.12p . 

Now we rearrange the columns of the determinant (|2.14p in the following way: 1st, 3d,. . . ,(2j — l)st, 
2nd, 4th,. . . , 2jth. To move the 3d column to the second place, we need one transposition. To move the 
5th column to the third place, we need two transpositions, and so on. At last, to move the (2j — l)st column 
to the jth place, we need (j — 1) transpositions. During all those transpositions the resulting determinant 

3- 1 i(i-i) 



will change its sign ^2 



times. Next we interchange the (j + l)st row with the 2jth row, the 



(j + 2)nd row with the (2j — l)st row, and so on. This will also require 



j(j - 1) 
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resulting determinant has the same sign as the initial determinant, so we have 
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Now from each (j + 2i — l)st row, i = 1,2,..., 



i + i 



we subtract rows £th, (£ + l)st, . . . , jth 



multiplied by to, ^2, ■ ■ ■ ,ti(j-i)-2i respectively. Then, from each (j + 2£)th row, i = 1,2, . 



subtract 



rows (£ + l)st, (£ + 2)nd, . . . , jth multiplied by t\, t%, ... ,t 2 (j-i)-xi respectively. As a result, we have a 
determinant, which can be represented in a block form: 



1 
A A 



fl 2j (h,g) 

where the upper triangular matrices Cq and Cj have the forms 
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Note that |C | = c£ > 0. 

In order to describe the matrices A and A, we establish a connection between the coefficients 6jS, CjS 
and t,;S. Multiplying (|2.ip by the denominator and equating coefficients, we get 



bk = Ck-jtj, k = 0,1,2, 



(2.17) 



i=0 

Taking into account those formulae, we obtain for the entries of the matrix A to have the form: 
fo if k < i + 1, 

— I — i — i. 

0>2i-X,k 

| °2(k-i) — 2_^ C 2(k-i-q) T 2q — 

••- 1 

where i = 1, 2, 



&2(fc— i) — c 2(fc— i-g)*2g = ^ C 2 ( fc _j_ g ) + i£2 9 -l if fc ^ £ + 1, 

9=0 

fi + i' 



7 7* "*1 



(2.18) 
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if fe < £ + 1, 
if ^ £ + 1, 



where £ = 1, 2, . 



From these formulae it follows that the matrix A can be represented as a product of two matrices: 
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Analogously, by (|2.17[) . for the entries of the matrix A, we have 
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where £ = 1,2,... 



Xhq = ^ C 2 (k-i-q)t2q+X if > £, 
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a-2i,k — < 
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a k < i, 

if k = i, 
if fc > i, 



where z = 1,2, 



Therefore, the matrix A can be represented as follows 
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It is clear that Aj is the j x j leading principal submatrix of the matrix %{R) defined in 
into account, we multiply both parts of the equality (|2 . 16[) by the following determinant 



E j 
C^ 1 



where Ej is the j x j identity matrix, to obtain 



c J n 2j (h,g) 



>0, 



A A,- 



or in the entry-wise from: 
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(2.21) 



The formulae (|2.18l) - (|2.2ip show that all the columns of the matrix A are linear combinations of the 
columns of the matrix Aj. Consequently, we are able to eliminate the entries in the lower left corner of 
the determinant (|2.2ip . To do this, from the 2nd column we subtract the (j + l)st column multiplied by 
Ci. Then from the 3 column we subtract the (j + l)st column multiplied by c\ and the (j + 2)nd column 
multiplied by C3, and so on. As a result, we obtain the determinant c^ 3 £l 2 j(h,g) to have the following 
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Since |Aj| = Aj(R), we get 
as required. 



n 2j {h,g)=<%A j {R), 



□ 



Now Lemma l2.4l and Theorem 12 . 2 1 imply the following criterion of Hurwitzness of real rational functions 
in terms of the coefficients of their numerator and denominator. 

Theorem 2.5. A real rational junction R defined in (|2.1[) - (|2.2I) is Hurwitz if and only if the inequalities 

n 2 (h, g ) > o, n 4 (h, g )>o,..., n 2n (h, g )>o 

hold. 

From (|2~TU1) and (|2~T3l) it also follows that 

fhn+2(h,g) = tt 2n +4(h,g) = Vt 2n +6(h,g) = ■ ■ ■ = 0. 

Lemma [2T4l and the formulae (|2.11j) imply the following corollary. 

Corollary 2.6. Let the polynomials h and g be defined in (|2.2| . Then the Hurwitz minors of the polynomial 
P(z) = (—l) m h(z)g(—z) satisfy the equalities: 



Aj(P)=n 2j (h,g), 

where n = deg P = deg h + deg g. 

Consider now the rational function 

F(z) = 



J = 1,2. 



(-1)" 



(2.22) 



R(- Z y 

where R is defined in (|2.1[) - (|2.2p . and n = r + m is order of R. 

It is clear that the function F is Hurwitz if and only if the function R is Hurwitz. This fact can be 
verified by the following relationship between the Hurwitz minors of the functions F and R. 

Theorem 2.7. Let the function R be defined in (|2.ip - (|2.2[) . and let the function F be defined in (|2.22l) . 

Then 

(2.23) 



A j (R)=t 2 j A j (F), j = 1,2, 



Proof. Let F(z) = , where f(z) = (— \) m g{— z), and q(z) — (—l) r h(—z). Consider the polynomial 
Q(z) = (—l) r q(—z)f(z), which is analogous to (|2.4I) . We have 

Q(z) = (-iyq(-z)f(z) = (-iy(-iyh(z)(-l) m g(-z) = (-l)™h{z)g{-z) = P(z). 
Now the formulae (|2 . 1 1 1) imply 

c^A^R) = A 3 {P) = A 3 {Q) = b 2 J Aj(F), j = 1, 2, . . . , 
bo 



that is exactly (12.23[) . since to 



co 



□ 
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Corollary 2.8. Let the polynomials h and g are defined in (|2.2[) . Then 

il 2j (h,g) = n 2j (f,q), j = 1,2,..., 

where f(z) = {-l) m g{-z), and q(z) = (-l) r h(-z). 

Finally, note that if one of the polynomials h and g is a constant, the determinants VL2j{h,g) become 
Hurwitz determinants (up to a constant factor) of the polynomial, which is not a constant. This follows 
from the formulas (12. lip . 
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